Universality in dissipative Landau-Zener transitions 
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We introduce a random variable approach to investigate the dynamics of a dissipative two-state 
system. Based on an exact functional integral description, our method reformulates the problem as 
that of the time evolution of a quantum state vector subject to a Hamiltonian containing random 
noise fields. This numerically exact, non-perturbative formalism is particularly well suited in the 
context of time-dependent Hamiltonians, both at zero and finite temperature. As an important 
example, we consider the renowned Landau-Zener problem in the presence of an Ohmic environment 
with a large cutoff frequency at finite temperature. We investigate the 'scaling' limit of the problem 
at intermediate times, where the decay of the upper spin state population is universal. Such a 
dissipative situation may be implemented using a cold-atom bosonic setup. 

PACS numbers: 03.65.Xp, 03.65.Yz, 33.80.Be, 74.50+r. 



I. INTRODUCTION 

A two-level system is never completely isolated 
resulting in dissipation, decoherence and entangle- 
ment [1]. Therefore, one primary task for experimental- 
ists is to manipulate and read out the internal state of 
the dissipative two-level system (qubit) with a high fi- 
delity. Often, this can be achieved by sweeping the two 
energy levels through an avoided crossing, a situation 
that occurs in a variety of physical areas such as molec- 
ular collisions [2] , chemical reaction dynamics [3] , molec- 
ular nanomagnets [4], quantum information and metrol- 
ogy [5-8]. For a constant crossing speed v this is known 
as the Landau-Zener problem [9-12] which can be solved 
exactly in the absence of dissipation. Naturally, it is im- 
portant to know the effect of the dissipative universe on 
the probability p(t) for the spin to remain in its initial 
state at time t [13-16]. Exact results [17, 18] are only 
available at zero temperature and in the limit t — > +oo, 
where the energy difference e of the two spin states is 
much larger than the bandwidth uj c of the environmen- 
tal bath. Typically however, uj c is much larger than 
the tunneling coupling between the two states A. Here, 
we rather focus on the experimentally relevant "scaling" 
regime at intermediate times, where the spin energies 
have not completely traversed the bath's energy band: 
A < e = vt < uj c with v > 0. To resolve the dissipa- 
tive spin dynamics, we develop a powerful numerically 
exact stochastic Schrodinger equation formalism (SSE). 
Compared to earlier SSE approaches [19-22], our method 
allows easier exact consideration of initial spin-bath cor- 
relations, which are crucial in the Landau-Zener context. 
It may also be applied to other many-body environments 
that can be represented in the form of a Coulomb gas 
such as the Kondo model. [23, 24] 

We prove that p(t) exhibits a universal decay in the 
intermediate (scaling) regime due to phonon assisted spin 
transitions. The size of the jump at the level crossing 
decreases for increasing dissipation and p(t) converges to 
the infinite time value only when t ~ uj c /v. We also 



derive an approximate analytical decay formula valid for 
slow sweeps at zero temperature, which agrees well with 
our numerical results. 



II. MODEL AND NOTATIONS 

Specifically, we study a two-level system coupled to 
a bath of harmonic oscillators (the spin-boson Hamilto- 
nian) [25, 26] 



— = —a 1 
h 2 



f 2 + yE^^ + M + E^- (1) 

k k 



Here, a x,z are the Pauli matrices, A is the bare tunnel- 
ing coupling and e the detuning. The bosonic oscillator 
operators have frequencies u>k and coupling constants 
We express the components of the reduced spin density 
matrix p(t) using functional integrals [25, 26] 

p(a f ,a' f ;t) = J Va(-) J Va'(-)A[a]A* [a']F[a, a'] , (2) 

where A[cr] is the amplitude for the spin to follow the 
path a(t) in the absence of the bath, and F[a, a'] is the 
real-time influence functional of the bath 
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written in terms of symmetric and antisymmetric 
spin paths ri(s) = ^[c(s) + c'(s)] and £(s) = 
|[cr(s) — °~'{ s )], respectively. The kernel func- 
tions Li(t) = Jp 00 duj </(w) sinwi and Li(t) = 
L°° dojJ(oj) cosojt cot\\ tvjj /2k bT are determined by the 
bath spectral function J(w) = hir Xf,S(ui — uju) and 
the temperature T . 

At time to — ^ — oo, the spin-bath interaction is first 
turned on, but the spin is held fixed in position <ii for 
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to < t < 0. The spin paths {a(t), a'(t)} in Eq. (2) 
are constrained to a(t) — a'(t) — a.- L for t < and to 
c(tf) — af,a'(tf) — a'j:. At t = 0, the bath is in the 
shifted canonical equilibrium state. For positive times, 
the spin jumps between the states {|t), II)} and the spin 
double path occurring in Eq. (2) can thus be regarded as a 
single path between the four states {|tt)> Itl), lit), III)}- 
If the path starts and ends in a diagonal ("sojourn") 
state { | ft ) > I II ) } an d makes 2n transitions at times 
t\ <t<2 < ... < i 2n along the way, it can be parametrized 

as £(t) = ET=i Si9(*-*i) and V(t) = ESo TM*-ti)- 
The variables {H x , . . . , H 2 „} = . . . , -£„} em- 

body the n off-diagonal ("blip") parts of the path be- 
tween the times t^m-i and t^m (m = 1, ...,n), and 
characterize the time spent by the path in the states 
{|tl>, lit)} such that {(t) = ±l,T](t) = 0. The variables 
{T , . . . , T 2 „} = {i] , -770, . . . , T] n } describe the n + 1 di- 
agonal (sojourn) parts in the time period (£2™, fem+i) 
during which rj(t) = ±l,£(t) = (here, we have m 
;</ and i 2Tl +i = i/). The path's boundary condi- 
tions then specify rj and r/ n . 

Inserting this general spin path ^{t),rj{t) into 
Eq. (3) and performing the time integrations yields 
Q1Q2 with 
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where Q12 are the second integrals of L 12 - The free 
spin-path amplitudes .A[cr].4* [ct'] give a factor to 
switch from a sojourn state rj to a blip state £ (and vice 
versa) as well as a bias-dependent phase factor H n = 
exp[i E?=i ^j s (tj)} witti s(t) = f*dt'e(t'). Altogether, 
the probability p(t) — p(|t), |t)j*) to nn d the system in 
state 1 1 ) at time t takes the form, 



K*) = i + E 



dt 2n • • • / ^ F n H n . 
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III. RANDOM VARIABLES 

We now proceed and decouple the terms bilinear in 
the blip and sojourn variables by Hubbard-Stratonovich 
transformations. Such a decoupling is useful since Eq. (6) 
has the Coulomb gas structure. [23] . Our formalism may 
thus be applied to other models which allow a Coulomb 
gas representation such as the Kondo model [24]. The 
resulting expression then suggests that p(t) can be ob- 
tained as a statistical average of a stochastic Schrodinger 
equation [19-23, 27]. 

For definiteness, we will now focus on the case 
of an Ohmic bath with spectral function J(uj) = 



77wexp(— uj/uj c ). It contains the viscosity coefficient r\ 
and a high-frequency cutoff uj Ci and we also introduce 
the dimensionless dissipation parameter a = rj/2irh. 
We like to emphasize that our method is able to solve 
for the system's dynamics at any temperature T. The 
bath correlation functions read Qi(t) — ?7tan -1 (u; c i) and 



tAtv sinh 



[25, 26] 



Q 2 (t) = %\n(l + LJ 2 c t 2 ) + ri\n 

In fact, to apply a Hubbard-Stratonovich transforma- 
tion to Eq. (5), we need to write Q2(tj—tk) in a factorized 
form Q 2 (tj - t k ) = 2[G + G ^ m fe)^(tfc)]. 
Since the kernel is translationally invariant, this can be 
achieved by a Fourier series expansion. To obtain only 
negative Fourier coefficients, we rather expand Qi(t) = 
Q 2 (r) - Q 2 (2) - a [ go + YZ=i 12 5m cos where we 

introduced the rescaled time r = t/t max , with i max be- 
ing the final time of our numerical simulation. Thus, the 
coefficients are Gq ~ " 



9o + ~<?2(2), G 2 k-i=G 2k =g k < 0, 
and the trigonometric functions read \E , 2fc-i=cos ^p, 
where k=l, . . . , m max /2. Decoupling the 
max Hubbard-Stratonovich transfor- 
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mations then results in 
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where the sum G=E™™' SVi. 
for 77i max — >oo, the integration over 
tributed Hubbard-Stratonovich variables reads J dS = 

TTl!!-i x f°° ^%s=e _s ™/ 2 , and we have introduced the real 

llm— 1 J— 00 T/27T 

function h{r) = J2Z=T SmV-aG m * m (r). 

We can proceed similarly with Q\ after separating it 
into a symmetric Qi(|f |) and an antisymmetric part Qi(t) 
in order to extend the sum to j < k. On the other hand, 
for zero detuning e = and a < 1/2 [25, 26], one can 
safely approximate Qi(t) ~ ry7r/2. This approximation 
becomes exact for A/ui c — > since the main contribution 
to the functional integral of Eq. (6) stems from spin flips 
with time separations larger than w" 1 . The finite bias 
case e ^ requires more consideration of the first sojourn 
as it accounts for the spin-bath preparation, which affects 
the long-time behavior ofp(t) [18] (see below). For e = 0, 
Eq. (6) reads 

p(r) = l + jdSY: ( , ) I dr 2n 

„ T2 n—l In 

x "'/ dTl E expect TyfcCfc+i] TT ex P[^ S 3^( T j)] 

(8) 

Without the summation over blip and sojourn vari- 
ables {£,j,r]j}, this expression has the form of a time- 
ordered exponential, averaged over the random variables 
{s,„}. This summation, however, can be incorporated 
into a product of matrices in the vector space of states 
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FIG. 1. (Color online): (a) P(t) as a function of t for various values of a, A = 1, ui c — 100, e = and T = 0. We checked that 
for a given a curves corresponding to different u c /A 3> 1 scale on top of each other in units of the renormalized tunneling rate 
A r = A(^-)° ! ' ( 1-a '. Quality factor of damped oscillations agrees with prediction n/7 = cot 2 {i- a ) f rom Refs. [25, 26, 28]. 

Results are obtained with m max = 3000, N = 5 ■ 10 4 . (b) P(t) for different temperatures T (here, h = ks — 1), dissipation 
strength a = 0.1, and other parameters as in (a). 
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e iira e ih(T) q q _ e -i-na e ih(r) 
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with A = ^(At m ^e-^^ 2 /^Q^+ G ^. Then, Eq. (8) 
becomes p(r) = J dS^^Te^ fo dsV ( s )|$;) which can be 
calculated by solving the stochastic Schrodinger equation 

i§ ? \Hr))=V(T)mr)), (10) 

with initial and final conditions | ) — (1, 0, 0, 0) T for 
N different realizations of the noise variables {s m }. Av- 
eraging the results gives p(r) = J2k=i ^[ k \ T )i where 
<I>i(t) is the first component of | $(t) ). Other com- 
ponents of the density matrix (2) can be obtained us- 
ing different initial and final conditions. In fact, the 
differential equations obey the additional symmetries 
lm$i = 0, $3 = $2 and $4 = 1 - $1, such that only 
three real- variables are independent. Since the evolu- 
tion is unitary (for e = 0) and + 2|$2| 2 = 1 is an 
integral of motion, we can introduce a classical unit- 
length spin S = (v / 2Re$2, V / 2Im ( E > 2, $1) that evolves 
according to dS/dr = H x S in a random magnetic 
field H = (cos h(r), sin h(r), 0), and from which we find 
p(r) = i(l + (S z (t))). Hence, the time-evolution of a 
dissipative quantum spin can be formulated as that of a 
classical spin in a random magnetic field. The quantum 
nature of the problem is hidden in the fact that spin ro- 
tations about different axes do not commute and through 
the averaging over random field configurations. 



IV. APPLICATIONS 

A. Spin dynamics at zero detuning 

To prove the feasibility of our method, we have com- 
puted the spin dynamics for zero detuning in the range 
< a < 1/2 for different temperatures T. We express 
T in units of A (hereafter we set h = ks = f). Results 
for P(t) — 2p(t) — 1 in Fig. 1 exhibit damped oscillations 
with the correct renormalized tunneling frequency of or- 
der A r = A(A)«/(i-«) for T < A r . The quality factor 
of the oscillations agrees with predictions from the Non 
Interacting Blip Approximation (NIBA) [25], field the- 
ory [28] and from the time-dependent numerical renor- 
malization group (TD-NRG) [30]. 

For intermediate values of alpha we are able to ac- 
cess the asymptotic long-time behavior of P(t), where 
|P(i)| <C 1, within our numerical approach. At T — 0, we 
find that the system exhibits exponentially damped co- 
herent oscillations as predicted in Ref. [28] and in agree- 
ment with recent TD-NRG calculations [30]. Note that 
the statistical error of our method scales like iV~ 1//2 , 
where N is the number of different realizations of the 
noise. We thus cannot resolve the existence of a small 
incoherent part at weak-coupling a In ^ <C 1 , which was 
predicted in Ref. 31 using a rigorous Born approximation. 

For increasing temperature, the coherence of oscilla- 
tions gets lost more rapidly, and finally for T 3> A,, 
we observe incoherent decay P(t) — exp[— t/r] with rate 

r_1 = £(l+l) *r[lB 2a t 25 ]' Note that our method 
gives reliable results over the full range of temperatures. 



B Dissipative Landau-Zener transition 
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FIG. 2. (Color online) (a) p(t) for a fast sweep with v/A 2 = 10, cj c /A = 200, a = {0.05,0.1} and T/A = {0,1,5} (here, 
h = ks — 1). We choose m max = 4000, N = 4 ■ 10 6 . (b) Slow sweep with v/A 2 = 0.5. Other parameters as in (a), (c) Fit of 
universal decay of p[e(i)] using Eq. (12) with a — 0.05 and single fit parameter C = 0.59. 



B. Dissipative Landau-Zener transition 

Next, we turn to the case of a Landau-Zener sweep of 
the detuning e(t) = vt with v > 0. We examine the sur- 
vival probability p(t) that the spin remains in its initial 
state if swept across the resonance. Neglecting the bath, 
this problem can be solved exactly [9-12] and one finds 
that pit) converges toward the celebrated Landau-Zener 
formula pi z = exp[— ttA 2 /2v) for t 3> A/v. 

A fundamental question is thus how this result is mod- 
ified in the presence of dissipation. Surprisingly, at zero 
temperature the bath does not affect the final transi- 
tion probability pi z in the limit t — > +oo, if the spin 
couples longitudinally to the reservoir via its a z compo- 
nent [17]. This limit, however, corresponds to very large 
times t 3> u> c /v where the separation of the spin ener- 
gies is larger than the bosonic bandwidth. In contrast, 
we explore the so-called scaling regime, where one first 
takes the limit uj c — > oo, holding A r t = y fixed, and only 
then considers y — ¥ oo. This limit is important because 
it allows the spin-boson model to exhibit universal be- 
havior [25, 26]. For large but finite lu c the scaling regime 
corresponds to an intermediate time regime where the 
spin energy separation e is smaller than w c but possibly 
much larger than A: A -C vt <C lo c . Phonon assisted 
spin transitions therefore still occur even though £> A, 
and the probability p(t) converges toward its final value 
pi z only for times of the order t ~ lo c /v. Note that this 
is in stark contrast to the non-dissipative (perfectly iso- 
lated) case where this convergence happens much faster 
for t ~ A/v. 

In the context of Landau-Zener transitions, the bath 
preparation affects the long-time result of p{t) [18]. Thus, 
it is important to consider the contribution of the initial 



sojourn exactly, as it accounts for the fact that the bath 
starts out from a shifted equilibrium state. It is given by 
the k = 0, 1 terms in Qi (Eq. (4)). We can incorporate 
this term by adding it to the height function 

- 2atan~ 1 [w c £ max (r - n)] . (11) 

Here, i max determines the time interval length of our 
simulation, and [0, r c ] ([r c , 1]) correspond to times be- 
fore (after) the level crossing. The fact that the height 
function now contains T\ forces us to explicitly per- 
form the Ti-integration in Eq. (8). We thus randomly 
pick a uniformly distributed t\ G [0,1], which deter- 
mines ft.(r, T\) as well as the initial state | $ Tl ) = 
-i(0,e i ^ ri ' Tl ),-e- l ' l(ri ' Tl \0) T . We then propagate this 
initial state in the interval [ri,l] according to Eq. (10) 
and calculate the survival probability as p(r) = 1 + 
($i(t)), where the average is over N choices of t\ and 
random variables {s m }. Here we set | <&(t < n) ) = 
in an individual run since ($i(t)) only accounts for the 
contribution of paths with at least one spin jump. In 
Eq. (10), the evolution is not unitary. 

In Fig. 2 (a), we check that p(t) converges toward pi z 
at long times t ^> w c /v for T = 0. For the large sweep- 
ing speed v/A 2 = 10 in Fig. 2 (a), we find that this also 
holds for T/A = {1,5}, since thermal effects only occur 
during the short period where < T [16]. In Fig. 2 (b), 
we show results for slow sweep velocity v/A 2 = 0.5, and 
observe clearly that the size of the jump in pit) at the 
crossing reduces with enhancing dissipation and temper- 
ature. Following the jump, we observe a decay of p(t) in 
the intermediate time regime up to t ~ lo c /v due to bath 
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mediated spin transitions. The final probability increases 
with temperature due to thermalization. 

We now derive an analytical formula describing the 
universal decay in the scaling regime at T = 0, which 
holds for slow sweeping speeds only. For large static de- 
tuning e 3> A r (but still e <C uj c ) the NIBA can be justi- 
fied [25] and predicts an overdamped exponential relax- 
ation with a decay rate T = ^p^y(e/A r ) 2a_1 . Inserting 

e(t) = vt and integrating dp/dt — —Tp(t) for a < 1/2 
yields 



p[e(t)} =Cexp 



-ttA? 



-(- e 



AaT (2a)v \A r 



(12) 



If we except the integration constant C, this formula con- 
tains only scaling variables, which shows that the de- 
cay is universal. It reduces to pi z in the limit a — > 
(with C = 1) and breaks down for times of the order 
t ~ oj c /v, where it becomes a function of the bare A 
again: p[e = uj c ] = Cexp[- 4a ^ a) J . In Fig. 2 (c), we 
show that the decay is indeed described by this formula 
for T = 0, v/A 2 = 0.5 and a = 0.05. Note that this de- 
cay does not occur at a = 1/2 [32] . We like to emphasize 
that our numerical method gives reliable results over the 
whole range of sweep velocities and temperatures. 



C. Realization with cold-atom quantum dot setup 

The intermediate (scaling) time regime A <C vt <C 
uj c might be accessed using the cold-atom geometry 
of Refs. [33-35]. It comprises a bosonic mixture of 
atoms in two hyperfine ground states a and 6, sub- 
ject to state-selective traps. One species forms a one- 
dimensional Bosc-Einstein Condensate (BEC), represent- 
ing the Ohmic reservoir, and the other species is trapped 
in a tight harmonic potential, operated in the collisional 



blockade limit, representing the "spin" . Coupling the dif- 
ferent species by Raman lasers, the system is described 
by Eq. (1) with A and e being proportional to the laser in- 
tensity and frequency, respectively. Using the parameters 
of Ref. [36], we estimate a = (— 1 + g a b/gaa) 2 ~ 0.06; 
K ~ \J Pa/gaa is the Luttinger parameter of the BEC, 
g a p = 2huj±a a p are the scattering amplitudes containing 
the transverse trapping frequency = 2tt x 67kHz and 
the scattering length a aa — 5.2nm. The value of a a t, must 
be tuned such that g a b *C g a a using optical Feshbach res- 
onances [37]. Choosing A w 100Hz and v s» lkHz/s, 
the intermediate time (scaling) regime occurs between 
0.1s<i<10s. 



V. CONCLUSIONS 

To summarize, we have developed a stochastic 
Schrodinger method to investigate the dissipative 
Landau-Zener problem in the scaling limit A/lo c <C 1 
at finite temperature. Assuming a < 1/2, we have shed 
light on an experimentally relevant intermediate time- 
regime where p(t) shows universal decay due to bath 
mediated spin transitions. Our results are relevant in 
quantum information, where fast quantum processes are 
more useful. Our method can also be extended to other 
many-body environments. 
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